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Abstract 

We study the triple differential distribution of i? ^ Xuii^, consistently including all 
perturbative and non-perturbative effects through 0(a^/3o) and 0(l/m^). The Fermi 
motion is parameterized in terms of a single light-cone function for each structure 
function and for any value of g^, accounting for all subleading effects. We discuss the 
problems and uncertainties related to the high-g^ tail and to Weak Annihilation effects. 
We work in the kinetic scheme, a framework characterized by a Wilsonian treatment 
with a hard cutoff /j, ~ IGeV. Our method is illustrated with the extraction of \ Vub\ from 
some of the latest experimental data, providing a detailed estimate of the theoretical 
uncertainty. 



1 Introduction 



The determination of the element \Vub\ of the CKM matrix in semileptonic B decays plays 
a central role in the search for flavour and CP violation beyond the Standard Model (SM). 
Currently, global fits to flavour violating observables predict [1] 

|Kfe|f/T/it = (3.44 ±0.16) X 10-3 

following the unitarity of the CKM matrix and assuming the validity of the SM. On the other 
hand, a direct determination of \Vub\ is possible based on the analysis of 6 — >• utu decays, 
either in the inclusive B — ^ Xulu or exclusive B —>■ Triu channels. Exclusive determinations 
rely on lattice QCD or light-cone sum rules for the corresponding transition form factors 
and have improved in the last few years. Two recent analyses give very consistent results: 

iKfelexci = (3.47 ± 0.29 ± 0.03) x 10-^ [2], iKblexci = (3.5 ± 0.4 ± 0.1) x 10"^ [3]. 

Inclusive decays based on a measurement of the total b ^ uiu decay rate potentially 
offer the most accurate way to determine \Vub\- They are described by a local Operator 
Product Expansion (OPE) in inverse powers of the b quark mass [1]. The OPE has proved 
quite successful in the analysis of the moments of various distributions in B ^ X^-lv, leading 
recently to the precise measurement of its dominant non-perturbative parameters, namely 
the b and c masses and the matrix elements of the relevant dimension 5 and 6 local operators, 
and to a 2% determination of | V^b| O [6]. 

In the case of charmless semileptonic decays experiments apply a series of cuts to isolate 
the charmless decays that tend to destroy the convergence of the local OPE. They introduce 
sensitivity to the effects of Fermi motion of the heavy quark inside the B meson, which are 
not suppressed by powers of 1/mb in the restricted kinematic regions. The Fermi motion is 
described in the OPE by a nonlocal distribution function, whose lowest integer moments are 
given by the expectation values of the same local operators we have encountered earlier. 

Fermi motion is of primary importance in another inclusive B decay, namely the radiative 
decay b ^ s'j. A dedicated OPE approach accounting for the relation to the nonperturbative 
i?-meson parameters extracted from B X^lv was developed and applied to the description 
of the photon energy moments with cuts in Ref. It proved quite successful in describing 
the available — *• + 7 data. The results of [Tj underline the importance of including 
subleading effects, going beyond the leading-twist description of Fermi-motion. Another 
advantage of the approach proposed in [7] is that it implements the Wilsonian version of 
the OPE with a 'hard' cutoff that separates the perturbative and non-perturbative effects 
[8] and reduces the significance of perturbative corrections. In this approach, sometimes 
referred to as the kinetic scheme, the non-perturbative parameters are also well-defined and 
perturbatively stable. The contributions of soft gluons are absorbed into the definition of 
the heavy quark parameters and of the distribution function. 

In this paper we develop an analogous approach for the case of the triple differential 
distribution in i? — ^> Xu decays. With respect to the radiative decays, there are however a 
number of complications due to the different kinematics. In semileptonic decays the invariant 
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mass of the leptonic system, g^, can vary up to M^. However the local OPE becomes 
problematic at high g^, where the effects of four-quark operators related to Weak Annihilation 
(WA) also show up. An accurate description of the physical spectra requires a careful 
inclusion of these effects. 

For what concerns perturbative corrections, the complete 0(a^/3o) corrections to the 
triple differential rate have recently been published [9\. To take advantage in our framework 
of this new results and of the well-known 0{as) corrections, we perform a new calculation 
of the real emission contributions with the Wilsonian cutoff both at 0{as) and 0(a^/5o). 
In this way, we include all perturbative corrections to the triple differential rate through 
O^alPo) in the kinetic scheme: to the best of our knowledge, ours is the most complete 
implementation of perturbative effects. The contribution of the 0(a^/3o) effects turns out 
to be numerically significant. In our treatment of perturbative corrections we do not resum 
Sudakov logs. Indeed, in the framework with the hard cutoff, soft divergences are absent 
by construction. The spectra still diverge at threshold due to coUinear divergences, but in a 
softer way. 

Another element of our approach that is common with Ref. [7] concerns the distribution 
function, namely the inclusion of its power corrections. We introduce the full finite-mb light- 
cone 6-quark distribution function. Power effects enter it through the power-suppressed terms 
in the OPE relations for its integer moments and we take them into account through order 
Aqcd Uni- At the level of power corrections the Fermi motion effects cease to be universal: 
they depend on the process, and for semileptonic decays they are function of and differ 
in the three relevant structure functions Wi-s{qo,q'^). We emphasize that we do not split 
the distribution functions into leading and subleading contributions. Dealing with the full 
finite-mb distribution functions we avoid calling upon a plethora of largely unconstrained 
subleading functions. The latter typically are increasingly singular in the end point, which 
is only an artifact of expanding in l/m?, rather than a physical effect. We study in great 
detail the dependence of the distribution functions on the assumed functional form. 

The significance of the effects of the Fermi motion proper in the differential distributions 
fades away at larger g^. However, at large g^ generic power corrections increase, and at some 
point even the integrated moments cannot be described by the OPE: for g^ approaching 
ml the decay process is no longer hard. This signals the emergence of WA effects. In our 
approach the change of the regime at certain g^ automatically manifests itself, at least as long 
as the 1 / ml effects are retained. To avoid the pathological behaviour of the OPE predictions 
in this kinematic region, we model the high-g^ tail in two different ways that preserve, for 
instance, the positivity of the differential rates. 

Our approach is implemented in a numerical C++ code and we illustrate it with the 
extraction of \Vub\ from some of the latest experimental results. Our results are compatible 
with the most recent Heavy Flavour Averaging Group (HFAG) average for \ Vub\ from inclusive 
decays jUj 

iKftlinci = (4.34 ± 0.16 ± 0.25) x 10"^ |Kb|mci = (4.31 ± 0.17 ± 0.35) x 10"^ 

where the two values refer to two different theoretical frameworks fTH [T3] . respectively, 
currently employed by HFAG. These values are a few standard deviations away from the 
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value preferred by the global fit to the unitarity triangle and by the exclusive determination. 
Our results, however, indicate slightly larger uncertainties, especially due to WA effects. 

The paper is organized as follows: in Section 2, after introducing some notation, we 
discuss the perturbative corrections to the triple differential semileptonic width in the kinetic 
scheme with a hard cutoff. In Section 3 we introduce the distribution functions and their 
convolution with the perturbative spectrum. In Section 4 we discuss different functional 
forms for the light-cone functions. Section 5 describes the problems encountered by the 
OPE in the high-g^ tail and presents two possible ways to handle that kinematic region. 
Our results are discussed in Section 6, where we extract \Vub\ from recent experimental 
data using the method developed in the previous sections. We also carefully discuss the 
various sources of theoretical uncertainty. After the Conclusions, the paper ends with three 
Appendices, containing certain details of our calculations. 



2 Perturbative corrections with a Wilsonian cutoff 

Our starting point is the triple differential distribution for B Xu^v in terms of the 
leptonic variables (go, (f', Ei) and of the three structure functions that are relevant in the 
case of massless lepton: 



d^T GlW, 



ub 



dq^ dqo dE^ Svr 
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2Ej - 2qoE, + ^ 



W2 + q\2Ee-qo)Ws] x 



xe[qo-Ee-^] 9{E,) 9{q') e{q. 



where go cind E^ are the total leptonic and the charged lepton energies in the B meson rest 
frame and g^ is the leptonic invariant mass. We will often use the normalized variables 



rrib ~ ml 



?o = — ' 1 = (2) 



where nif, is the b quark mass. Within its range of validity, the local OPE allows us to 
separate perturbative and power suppressed non-perturbative contributions to the structure 
functions 

ly,(go,g^/i) = rubif^r [Wr\qo,q^f^) + W^^l^o,g^/i)] (3) 

with ni^2 = — 1 and = —2. The quantities W[°^ contain the power corrections of the local 
OPE: their expressions through 0{l/ml) are quoted in the Appendix [Bl In the context of 
the OPE with a Wilson cutoff /x, the separation between Wf^^^ and Wf°^ is controlled by /i 
and both contributions are /^-dependent. The /^-dependence, of course, cancels out at each 
perturbative order in inclusive quantities like the go-moments of Wi(go, g^, /x). 

As already mentioned in the Introduction, we absorb the contributions of soft gluons 
in the definition of heavy quark parameters and of the distribution functions. Physical 
quantities are in principle independent of the cutoff. The presence of the cutoff introduces 
several modifications in the perturbative structure functions. They have been studied at 
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< (1 - 2 r])' 





Figure 1: The real gluon emission energy in the presence of a cutoff, with t = 2Eg/mb. In the 
on-shell decay amplitude t is constrained to lie within the black solid line. The red horizontal line 
represents the cutoff t > 27]. The left (right) plot is at low (high) g^. See the text. 

= for radiative decays in [7] . In our case the structure functions take the form, through 

oia%y. 



Vi [Q ,V) +Cf^^V^ '{q ,V) 



6{l + f-2qo) 



+ Q 



71 



Bl '{q ,V) + Bl '{q ,r]) 



e{l + (f- 2go) 



TT 



5'(l + g'-2go 



(4) 



where 



T] = fJ'/mb. 



(5) 



We always assume < 1] < 1/2, and i3q = 11 — |nj with Uf the number of hght active 
flavours. In the numerics, we will set = 3. The derivative of the Dirac's delta in (jl]) is 
taken wrt to its argument. 

The normalization of the Wi is such that 



(6) 



We have computed the real gluon emission contributions i?- ' {qo,(i'^,v) restricting the 
phase space integration in b —>■ u W* g to gluons with energies larger than the cutoff fi. The 
effect of the cutoff is to remove the infrared divergence, softening the divergence of the form 
factors at the endpoint, where collinear divergences are still present. The calculation of BLM 
corrections has been performed using the technique with a massive gluon and integrating over 
the gluon mass, a standard procedure described, for instance, in [H], that can be extended 
to compute and resum 0(a"/?o~^) corrections. We have also reproduced all results of the 
analogous calculation performed in [7] for the case of 6 ^ 57. 
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Figure 2: Different domains for R\ ' in the {qo,q^) plane. The solid and dashed straight lines 
represent (1 + g^)/2 and w{(f,ri), respectively. Regions I, II, and III comprised between these 
two lines and the lower bound are shown in dark gray, light gray and red (see text). We 
have used /x = 1.2 GeV, = 4.6 GeV. 



The functional form of the real contributions Rf^\qQ, (p, if) depends on the region of the 
parameter space: 

Rf\%. e, V) = Rf\%, f) e {w - go) (l - 2r^ - v^) + ^'^('?o, v) {% - w) (7) 
where 

^"^^^''^) = ^-^+ 2(l-2r^) 

and R^''\qo , g^) are the real emission contributions in the on-shell scheme (without cutoff, 
/i = 0) calculated at 0{as) [15] and 0(a^/3o) [9J. Eq.(l7]) can be understood from Figs. [H 
where the integration range for the gluon energy Eg in the b rest frame is represented in 
terms of t = 2Eg/mh and go at fixed (f. The on-shell kinematics implies 1 — go ~ \/ qi~ q^ < 
t < 1 — go + Y^go — g^. At low g^ (left plot) the cutoff is irrelevant for go < w{q'^,ri) while 
it modifies the spectrum for w^q"^,!]) < % < (1 + g^)/2. For large enough values of g^, 
g^ > (1 — 2//)^, the energy of the gluons is always below the cutoff unless w^q"^,!]) < % 
and the spectrum is affected at all go (right plot). At even higher lepton invariant mass, 
g^ > (1 — 2r]), gluon emission is completely inhibited. One therefore identifies three regions 
in the (go, g^) plane that are displayed with different colors in Fig. [2j (I) where the cutoff 
does not modify real emission (dark gray), (II) where the cutoff modifies the real emission 
(light gray), (III) where the presence of the cutoff inhibits completely real gluon emission 
(red). 
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Figure 3: The real gluon contribution i?^ ' (dashed curves) and the sum n\ ' + asPo/nRi 
(solid curves) as functions of go at (p = 0.2 with (red) and without (black) a cutoff ji = 1 Ge\/ 
on the energy of the gluons. 



Explicit expressions for R^^\qo,q'^,ri) are reported in the Appendix A. We do not give 
those for Rf'\qo,q'^,i]) which are relatively lengthy, but they can be found in the computer 
code. Fig. [3] shows the go dependence of Wi with and without the cutoff. Here and in 
all numerical examples of this section, we assume /x = 1 GeV, mb(l GeV) =4.6 GeV. As 
anticipated, the structure function is unchanged for go < w, and diverges less severely close 
to the endpoint. Therefore we do not perform any resummation, unlike [121 [131 EH]- 

A direct calculation of the virtual contributions V- ' in the presence of the cutoff fi 
is more cumbersome. However, since their expressions in the absence of the cutoff (on- 
shell scheme) are known from |15l [9], one can infer their expressions at arbitrary fi from 
the requirement that physical quantities be independent of the cutoff at each perturbative 
order when both power and perturbative corrections are consistently included. In particular, 
following the argument given in the previous section, we require the /i-independence of the 
integral over go of each structure function: 

dqoW^iqo,q^O)= dqoW,{qo,q\ ^) + 0{al). (9) 

-OO J — OO 

In order to extract the /i-dependence of l^, we therefore consider the various sources of 
cutoff dependence in Eq. The renormalization of the non-perturbative parameters and 
of mf, is known to induce a cutoff-dependence of the form [S], [T7] : 

/X^(0) = - [/i^(/i)]pert (11) 

PD(0)=pL(/^)-[PD(/")]pert, (12) 
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[^(A^)jpert 
[/^'(/^)]pert 



4 as{mb) 

3 TT 

Q;s(mfc) 2 

Op fJ' 

TT 



1 + 
1 + 



TC 



1 + 



271 

2tt 



irib H 

m — 

2/i 3 

m;, 13 

m 

2/i 6 



2/i 



(13) 



It is worth recalling that the perturbative shifts in Eqs. (11011121) are not only conceptually, 
but also numerically quite important: using the current best experimental determination 
/i^(lGeV) ~ 0.40 GeV^, for instance, we see that the perturbative shift [/i^(l GeV)]pert 
amounts to almost 40% of that. We now write the perturbative contributions to the structure 
functions in the following way: 



^1,2) 



(1,2), 



(14) 
(15) 



where V^^^'^'*(^^) are the soft- virtual contributions in the on-shell scheme (without cutoff) 
calculated at 0{as) [15j and O^a^Po) [2] (see Appendix A). 
The matching condition (Q is satisfied if and only if 
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dq^5Rf\qo,q\v) 9{l + f-2qo 



(16) 



where the range of the integral on the rhs is fixed by the decay kinematics. It is best 
performed numerically, while the S'l^\(f',r]) are given by 



Cpasimb) 



TT 



IT 



m 



2 



[A(/i)] 
d 



pert 



+ 



pert 



2mb 
_d_ 

pert gp3^^ 



(0),tree 



dmh 



M 



{O),pow 



if) ■ (17) 



I pert Qp2 

The tree-level and power corrections contributions to the zeroth moments, Mj^'^^'^^'^'^ and 
, can be found in Appendix B. The resulting factors S^''^ read: 

,2 TQ^^2_2) 



(1) 



8 .2 ^ 2 

3? V + ^—V , 



(1) 



0, 



5f : 



32g^ , 13(g^-2) , ^(i) ln2r/ 
9 ' 36 ' ^ 2 ' 

0, 



(1) 



•3^+9^' 



S. 



(2) 



32 2 , 
T^+9^^ 



s. 



(1) In 2r] 



^ These expressions actually refer to the asymptotic value of /x^, namely in the infinite limit. In 
general we employ a definition of and of the other OPE parameters at finite m{,, but we actually choose 
to neglect 0{iJ,^) terms in [/x^(/z)]port and [/Zg.(^)]pert in order to be consistent with Ref. [18] and with the 
way the global fit for the determination of these parameters is performed in [B] . This amounts to an ad-hoc 
perturbative redefinition of /i^ and /i^. 
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Figure 4: Virtual contributions to V^^\q'^,fi) (left) and V^'^\q'^,fi) (right) and their pole mass 
scheme counterparts V^^'^^ = Vl'^'^\q'^ ,0) (dashed lines) as functions of for /i = 1 GeV and 
tts(^fe) = 0.22. The blue, red, and magenta lines correspond to i = 1,2,3, respectively. 

The complete virtual contributions to the three structure functions are shown in Fig. IHfor 
fi = 1 GeV, at 0{as) and 0{alPo). There is a strong suppression of the virtual contributions 
with respect to the case without cutoff and the BLM corrections are typically smaller in 
relative terms, as it can be expected since the cutoff increases the typical gluon energy. It 
is worth noting that, following our discussion of Figs. there are three regions of in 
the calculation of the integrals of 5Rl ' in f|T6|) . This is illustrated in Fig. [5] for the case of 
the virtual contributions Vi. For > (1 — 21])'^ terms non-analytic in yU = appear: the 
radius of convergence of the /^-expansion decreases for increasing and one does not expect 
the //-independence to hold when one includes only a few higher dimensional operators, as 
done above. This is actually related to the poor convergence of the OPE at high g^ that 
we will discuss more in detail later on. In practice, our choice of a cutoff close to 1 GeV 
implies a rather small value [1 — 21])'^ ^ 0.32 for the first threshold in g^. A high cutoff might 
indeed magnify at moderate g^ the contribution of higher dimensional operators that are 
inevitably important in the high-g^ tail. For this reason in our implementation we keep the 
complete ?7-dependence of SVi, but one can examine different options to study the theoretical 
uncertainty. 

There is one point left to clarify in Eq. (jlj), namely the presence of a 6' term. This is 
related to the difference that occurs between the kinetic mass of the b quark and the rest- 
energy that determines the kinematic end-point in the decay of the heavy quark |7J. This 
difference implies a finite perturbative shift of the endpoint, which manifests itself in Eq. (jl]) 

(1 2) 

as a derivative of the Dirac delta. The values of ' can be calculated explicitly, but again 
we choose to infer their expressions from the //-independence. To this end we require the 
cutoff-independence of the first go moment (the only moment affected by the presence of a 



8 



Figure 5: Virtual contributions to the structure functions at 0(«s), Vi^^\q'^,^) (thin lines) 
and through 0{all3o), v}^^ + ^vl"^^ (thick lines) as functions of for = 1 GeV and 
as(mfe) = 0.22. The blue, red, and magenta lines correspond to z = 1,2,3, respectively. 



S'). A procedure analogous to the one described for the virtual corrections yields: 
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?(2)/;.2 
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2+ni 



pert 
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pert 



2mh 



m 



2+ni 



d 



d 



r{l):POw/;52 



Cpas 



I port ' iru^t-Jlpert Qp3^^ 



^(0),tree^^2) 

(19) 



dqo % 



TT 



where the power corrections to the first central moments of the structure functions, 

are given in the Appendix B. In the region < (1 — 2?])^, expanding in rj up to 0(r7^), we 

find 



(1,2) 



B. 



(1) 



B 



(1) 



^(l-g2) ,^ + 4r/2 + l(l + 5g2) r/^, 
^(l-g2),^ + 2r72^^(-7 + 5g2) r/l 



(20) 
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Figure 6: The coefficients Bi ' as functions of (f for /i = 1 GeV. The blue, red, and magenta 
thick lines correspond to i = 1,2,3, respectively. The thin black lines are the corresponding r] 
expansions given in Eq. (I2U]) . 



The Bi are given as high precision approximate formulas that reproduce the numerical 
results within 1%. 

As noted above, non-analyticity in rj is significant at high (f . This is again illustrated in 
Fig. [6] for the case of Bf^. The figure compares the full ?7-dependence with the expansion 
in 7] through 0{rf). For > (1 — 27])"^ there is no reason to expect Bi to be well described 
by an expansion in rj. It turns out that non-analyticity is a minor effect for i = 1, 3 up 
to (f' ~ 0.8, while for z = 2 it becomes relevant already at (f ~ 0.5. On the other hand, 
W2 is kinematically suppressed in the triple differential rate at high q^. In our practical 
implementation of perturbative corrections, we employ only the expanded formulas fl20l) 
for Bl . the consequent mismatch at high q manifests itself as a mild //-dependence of 
the physical quantities that we compute and will be taken into account in our estimate of 
theoretical uncertainties. 

3 Distribution functions in B ^ Xuli^ 

The resummation of the leading twist effects in semileptonic B decays close to threshold has 
been first studied long ago [19]. The structure functions are expressed as the convolution 




1.444(1 - f)7] + (4.691 - 0.0544g2 + Qm2q^)7f 
-(0.7517 - 4.353g2 + 5.612g^)r/^ - B^^^ 
0.7219(1 - q^)r] + (2.353 - 0.044g2 + {).{]hq^)rf 




-(2.772 - 2.65g^ + 2.Sf)rf - B^^ 
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of their tree-level expressions with the light-cone distribution function F{k^) whose support 
lies below A = Mb — rrih, where Mb is the B meson mass: 



m. 



dk+ F{k+) Wl 



(0) 



90 - — 1 1 - — 



mf 



(21) 



where 



k+ is the light-cone component of the residual h quark momentum, and F{k+) is the dis- 
tribution function. Eq. (12T1) is valid at the leading order in l/rrib and does not include 
perturbative contributions. 

The main properties of the distribution function F{k^) in Eq. fl2T]) are well-known: while 
it cannot be presently computed from first principles, its moments follow from the local OPE 
because they are related to the go-moni^nts of the structure functions Wi. The distribution 
function is universal, i.e. independent of the considered structure function and shared by 
inclusive radiative and semileptonic B decays. Moreover, since the leading-order moments 
of F{k^) are independent of g^, the distribution function does not depend on g^. 

Perturbative corrections can be included in the leading twist formula Eq. (pTj) . by using 
in the convolution the short- distance perturbative structure functions that contain gluon 
bremsstrahlung and virtual corrections, instead of the tree-level W-^\qQ,q'^). 

The phenomenological and conceptual importance of subleading contributions to Eq. ( 1211) 
has been repeatedly stressed in the last few years \T, '20^. In order to proceed beyond the 
leading order, we first of all modify Eq. fl2Tl) into 



Wi{qo,q' 



m 



dk+ F{k+) Wl 



pert 



^° 2 V rribMB 



(22) 



The advantage of the latter representation is that it automatically yields the correct hadronic 
endpoint for go arbitrary g^: since k^""^ = A, the maximum value for go is 



q'o 



rUbM, 



B 



+ 



mt 



2mb 



Ml + q^ 
2 Mb 



(23) 



Following [7] we do not split the distribution function into separate leading and subleading 
components. Rather, we define the finite-m;, distribution function assuming the form of the 
convolution in Eq. (1221) to hold at any order. 

Clearly, such a distribution function is no longer universal: there must be a distribution 
function for each structure function and none of them corresponds exactly to the one describ- 
ing the radiative decays. Similarly, since their moments including power corrections depend 
explicitly on g^, the distribution functions depend on g^. Finally, like the OPE parameters 
and the perturbative corrections, the distribution functions depend on the infrared cutoff 
/i: by construction, they are designed to absorb all infrared physics characterized by energy 
scales below /i. The form (l22l) of the convolution amounts to their perturbative definition. 
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We will therefore work with three distribution functions Fi[k^, g^, fi) that, after including 
the perturbative corrections according to Eq. fl22]) . must lead to /x- independent distributions. 
Thus, our generalized convolution reads 

^ ^ ^ (24) 



W.{qo,q')=m^^{fi) J dk+ F,{k^,q^fi) 



1 



q^ 



rribMj 



B 



As mentioned already, all the available information on the distribution functions is en- 
coded in their momentqj. 

They can be extracted by matching with the OPE predictions for the go-^ioKieiits of the 
structure functions, known through 0(a^/3o) and Since the perturbative corrections 

to the Wilson coefficients of the power suppressed operators are not known, the moments 
of the Fj(/c+, g^, /i) can be determined neglecting perturbative corrections. We now consider 
the following moments of the structure functions 

j dqoiqo - ar W,{qo, q') = m^' J dk+ F,{k+, q') j dqo{qo - af ^^^(go - fc+ y, g'), (25) 



where we have left all /i-dependence implicit. We have also employed 

a = "^^ ' , A = 1 



2 mfe 



rribMB' 



(26) 



and replaced W[^^^ with its lowest order term W-^\ The Ihs of Eq. ( !25l) is calculated including 
power corrections and reads 



dqoiqo -ar Wt\qo,q') + Wriqo,q 



m 



Ui+n+l ( j(n),tree 



+ 1, 



{n),pow 



with the /,^"''*'s given in the Appendix B, while the rhs becomes 



m 



rti+l 



r(0),tree 



J dk+ klF,{k+,q'), 



(27) 



(2^ 



where we have used the vanishing of all central moments of the tree-level structure functions 
is also given in App. B. The moments of the distribution functions are 



(0),tree 



for n > 0. 
therefore given by 



dk+ klF,{k+,q' 



2 rrih 



r{n),pow 

i 

T(0),tree 



(29) 



The functions Fiik+^q^) are normalized to 1, up to small power corrections. The q^ depen- 
dence of the first two moments is shown in Figs. [3 The importance of subleading contribu- 



^The photon spectrum of i? ^ Xg^ is currently measured with a good accuracy at the B factories for 
> 2.0 GeV. These data are important for the precise determination of nib and of the OPE parameters. 
However, the underlying distribution functions in semileptonic and radiative decays differ at the subleading 
level, and some of the subdominant effects do not seem to be under control [H]. Therefore, we believe 
that the photon spectrum in i3 ^ Xg^ should not be used directly to determine the distribution function 
of semileptonic decays. On the other hand, as illustrated in [7] and its phenomenological applications, the 
photon spectrum can be accurately predicted, using our method, for not too high. 
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Figure 7: Dependence of the first and second moments of the distribution functions Fi{k^, g^, /x) 
on g^. The blue, red, and magenta curves correspond to i = 1,2,3, respectively. 

tions is apparent: the first moments for the three structure functions are quite apart from 
each other and we observe a strong g^-dependence at high g^, where they even diverge. The 
second moments of Fi^^ decrease with increasing g^ and have a zero at g^ ~ 0.7 — 0.8. The 
variance of these two distribution functions decreases for growing g^, until it reaches negative 
values and the very concept of (positive definite) distribution function becomes meaningless. 
This is not surprising: as already noted in ^19j, the light-cone distribution function of the 
heavy quark cannot describe the semileptonic decay at high g^. We will come back to this 
subject in a dedicated section. 

It is possible to use alternative forms of the convolution of the short and long-distance 
contributions, which differ at the subleading level from Eq. ( 12^ . For instance, one could use 
ml or instead of mbMs in the argument of W^^^^. However, this would simultaneously 
redefine the power corrections to the moments of the light-cone function, and the consequent 
change in the distribution functions would largely compensate the change in the convolution 
at the level of observable structure functions. Therefore, adopting our procedure the choice 
of convolution is not an ad-hoc assumption. 

Once the distribution functions are required to respect the OPE relations for their mo- 
ments with power accuracy, the only element of arbitrariness concerns the ansatz that is 
employed for their functional form (see next Section). The latter is not specific to the treat- 
ment of the subleading-twist effects and is present already in the leading-twist analysis. 
Varying the functional form will allow us to estimate the associated theoretical uncertainty. 
The convolution ( 12^ can now be used to compute the structure functions. It is useful to 
make the ^-function contained in F{k+) explicit: 



Fi{k+, g^ /i) = Giik+, q\ /i) e{K - k+), 



(30) 
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Figure 8: The convoluted W3 at = (left) and = 10 GeV^ (right) for /i = 1 GeV. The 
dashed, dotted, and solid curves correspond to lowest order, next-to-leading order, and 0{alPo) 
in the perturbative corrections, respectively. 



and write the detailed form of the convolution as 



m 



1 

A 



2go - 1 - ,2 



+ 



1 CpOs 



A2 



TX 



B] '{q ,7]) + B:> '[q ,rj) 



TT 



71 



A 



dK Gi{K,q^,n) 



G 



A 



, /2go -l-f 



A 



/ / TT 



A 



A 



go 



(31) 
-K,q' 



Here we have employed A = A/mfe. In Fig. [8] we show the go-dependence of the convoluted W^, 
at = and g^ = 10 GeV^ for /i = 1 GeV for the exponential ansatz discussed in the next 
section. The plots compare the structure function calculated at the leading order (dashed 
curves) and at the next-to- leading order (dotted curves) in a^, and including also the Oi^alPo) 
corrections (solid curves). We observe that the width of the convoluted structure function 
shrinks significantly between g^ = and g^ = 10 GeV^. The main effect of perturbative 
corrections at small g^ is a broadening of the shape, due to real gluon emission. At higher 
g^, instead, real emission is progressively inhibited. In the example considered, g^ = 10 GeV^ 
lies just below (1 — 2?]) (see Fig. [2]). Therefore the only appreciable perturbative effect 
is a shift of the peak somewhat towards higher go values, and is driven by the 5' term in 
Eq. (jlj). The Wi^2 structure functions behave in a very similar way. 



4 Functional forms 

We have seen that only the first few moments of the distribution functions are known. 
They are given in terms of matrix elements of local operators that are measured in 6 — c 
semileptonic and radiative decays. The modelling of QCD dynamics in the threshold region 
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Figure 9: A comparison of the basic functional forms in Eqs. ( J32H36I) for the distribution function 
g^, /i) after using the first two moments at = and = 8GeV^. 

therefore requires an ansatz for the functional form, which must comply with the constraints 
coming from the first few moments. Several functional forms have been proposed in the 
literature. They have to be exponentially suppressed at large negative fc+ and to vanish at 
the endpoint /c+ = A. We will further assume the positivity of As suggested by its 

probabilistic interpretation, the primordial Fi{k^) should be positive definite, but radiative 
corrections can potentially change its sign, in analogy with the parton distribution functions 
of deep inelastic scattering. In our approach with the Wilson cutoff effects of this kind are 
excluded in the negative tail and we choose to neglect them altogether, assuming positive 
distribution functions. 

The basic two-parameter functional forms proposed in the literature [22l [T3t [7] are 

F{k+) = N{A- k+ye^''+ e{A- k+) (exponential) (32) 

F{k+) = N{k-k+fe-^^~^-''+^^e{K-k+) (gaussian) (33) 

F{k^) = N ^^^^^"1^"/^^°^ .1 ^(A - K) (hyperbolic) (34) 



, 2 

-ai A— A;j_ + 



F{k+) = Ne y ^ ^-'+/ e{A-k+) (roman) (35) 

The parameters a,b > are fixed by the first two normalized moments of F{k^) and can be 
easily found numerically. A comparison of the resulting distribution functions Fi(fc_|_, g^, /x) 
at = and g^ = 8 GeV^ is shown in Fig. [91 The numerical inputs for the moments 
are taken from [6j. The results for the other two light-cone functions are very similar. The 
hyperbolic and exponential forms are almost indistinguishable in the plot. Indeed, once the 
first two moments are fixed, the shape of the various curves is determined to a large extent 
by the properties of the tails. In this respect, there is not a large variety in ( I32i) -( l35l) . It is 
not difficult to find alternatives: for instance, one could modify the exponential form as in 

F{k+) = N{A- k+Y e*(^-'=+)^ e{A - k+), (36) 
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Figure 10: Exponential form for Fi(A;+, g^, /i) retaining only the leading contributions to its 
moments (black curve) and including all known subleading effects at = (red) and at = 
lOGeV^ (magenta). 

which is also displayed in Fig. [9] (dashed curve) and is characterized by a higher negative tail. 
It is interesting to see the effect of subleading corrections to the moments of the distribution 
functions: this is illustrated in Fig. [TOl where Fi{k+, g^, /i) is computed using the exponential 
ansatz retaining only the leading contributions to its moments, and including all the known 
power corrections. 

The need for more flexibility in the choice of the functional form can also be understood 
by comparing the third normalized moments of the curves displayed in Fig. [91 For a given 
value of g^, they are remarkably close to each other. Let's consider g^ = 0: while the gaussian 
form yields —0.054 for the third normalized moment, exponential, hyperbolic, and roman 
forms cluster around —0.071, and that in yields —0.079. On the other hand, the OPE 
predicts the third moment to be — p|)/3, up to unknown 0{l/mh) and 0{as) corrections. 
Using the results of the fit in [6], this amounts to —0.058 ± 0.007, in rough agreement with 
the values mentioned above. The situation does not change significantly with g^, but the 
spread between different forms decreases further. 

Clearly, a mismatch between the third moment of the two-parameter forms we have 
considered and the OPE prediction might signal either large and possibly g^-dependent 
subleading contributions to the OPE, or that the specific functional forms are disfavored 
by present data on inclusive semileptonic moments. A line is difficult to draw, but none 
of the basic forms in fl32l - [36|) appears particularly disfavored. Rather than using the third 
moment to constrain the distribution functions, we have used it to gauge the diversity in the 
functional forms we employ. 

A straightforward generalization of the above basic forms consists in multiplying them 
by a distortion factor D{k+) that can be a polynomial, positive in the (— oo. A) domain, or 
a more general positive function. In principle the distortion can depend on g^. Of course 
not all the distortions D{k+) yield acceptable solutions. In particular, we discard shapes 
with two or more maxima, but keep those with an inflection point. We are particularly 
interested in distortion factors that modify the third moment. This can happen when the 
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Figure 11: Left: a sample of possible ways to distort the exponential form (red curve) in Eq. (1320 
compatible with the first two moments for the distribution function Fi{k^,0, IGeV). Right: the 
same but considering other basic functional forms. 

distortion enhances or suppresses the lower (negative k+) tail, or, to a lesser extent, with an 
asymmetric and bounded function, like a constant with a sinusoidal perturbation. In general 
it is difficult to find distortions that modify significantly the third moment. A representative 
but not exhaustive sample of distribution functions Fi{k+, = 0), based on the exponential 
form and complying with the first and second moment constraints, is shown in Fig. [TTl All 
the functions in the sample satisfy the first two moment constraints and yield third moments 
that differ by less than 30% from the basic exponential form. Since there is no reason of 
principle to discard any of them, they will be used, together with similar ones, in our study 
of the theoretical uncertainty. Changing also the basic functional form one gets the plot on 
the rhs of Fig. [TTl with almost a hundred different forms. 



We have seen that the formalism developed in Sec. 3 cannot be applied at high g^. The 
relevance of Fermi motion subsides at high and Fig. [7] shows that the variance of the 
distribution functions F13 becomes negative at ~ 0.7-0.8, as a consequence of O(A^) 
effects. While the leading contributions to the n-th central moment of the structure function 
Wi, j^^^ / j^'^^^^'^^ ^ is suppressed by (1 — g^)" at large g^, subleading contributions in l/rrii, have 
a weaker suppression that enhances their weight at high g^. 

We also recall that our discussion of perturbative corrections points to larger uncertainties 
in the high g^ region. We have seen in Sec. 2 that the expressions for the perturbative 
corrections become non-analytic in rj = n/nib for g^ > (1 — 2?7)^. This is a manifestation 
of the same problem mentioned above: at high g^ the contribution of higher dimensional 
operators is no longer suppressed. 

A simple way to visualize the problem is to plot the width of the physical range in go as 



5 The high-(?^ tail 
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Figure 12: The width H of the go range expressed in GeV as a function of for = 4.6 GeV. 
The gray band corresponds to H < a/ /U^/3(1 — g^) (see text). 

a function of g^, as in Fig. [121 The OPE structure functions are smeared over a region in go 
of width S = {m}j — /2mi,. H shrinks rapidly with g^, reaching a size comparable to 

^QCD for ~ lOGeV^. The range S can also be compared with the typical width of the 

I 2 

structure function, which at the leading order in l/m;, is y ^(1 — (f) ~ 0.36 (1 — (p) GeV. 

As shown in Fig. [121 they are equal for g^ ~ 11 GeV^ and for higher g^ the width of the 
physical go-range is smaller than the width of the structure functions, a situation that makes 
all local OPE results unreliable. 

It is indeed well-known that for very large g^ the local OPE fails to provide a reliable 
description of the semileptonic decay, as the dynamics in that region is not characterized by 
a short- distance scale [191 [23]. The Wilson coefficients in the power expansion, for instance, 
grow like 1/(1 — g^)". The origin of this enhancement at high g^ can be understood by 
calculating the g^-differential rate as the integral over go of the double differential rate. 



dV 



dqo dcp- 



ql-cP[^^W, + {ct,-(f)W,]. (37) 



It is related to the presence of the square root, a non-analytic term that follows from the 
kinematic constraints on the electron energy — the first 6'-function in Eq. Expanding 
the square root in go around the endpoint go = (1 + g^)/2 one gets 

l-g2 - (_2)«6„(g2) / l + q 



q5-Q' = — — Q _ ^2)2n-i yio — —j . (3^ 



52 _ ^2 = ' 

2 

n=l 



Clearly, the convergence radius of the series gets smaller at higher g^ and, in particular, 
higher orders become more and more important. Replacing the square root in (l37j) with 
its expansion (l38l) and integrating over go one sees that non-perturbative contributions are 
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organized in the series 



nib 



(39) 



dq^ (1 _ g2)n-2 

n=l ^ ' 

and one can verify that unknown 0(l/m^) corrections become comparable to 0(l/m^) cor- 
rections at ~ 0.7. 

Since the contributions of higher dimensional operators become more and more singular 
at high g^, we see from f l39p that their contribution to the h^qlv total width is singular for 
n > 3. In particular, the Darwin operator contribution to the q^ differential rate is {<f ^ ^) 



dq^ 6m 



96 

20 f + 66 + 48 q^ + 74 



1-q' 



which upon integration over leads to a logarithmic singularity that can be regulated by the 
quark mass flU{ [23] . This unphysical singularity is removed by a one- loop penguin dia- 
gram that mixes at 0{as) the Weak Annihilation four-quark operator = —Abl'yul u^'yb^ 
into the Darwin operator. The cancellation is discussed in detail in [23]. Defining 



B^^ = {B\01JB), = 327rVm3, (41) 

one finds that the sum of the contributions of the WA and Darwin operators to the semilep- 
tonic total width is 



6T 



CwA 5wA(/iwA) - ( 8 In ^ - 5 ) ^ + (^(« 



(42) 



which is free of singularities. The scale /Xwa is the MS renormalization scale of the WA oper- 
ator. The constant accompanying the logarithm in Eq. (H2l) depends on the renormalization 
scheme; we have employed here the same scheme as in Ref. [25]. At the level of the 
differential spectrum of Eq. fj^Oj) . the singularity 1/(1 — g^) is replaced by 1/(1 — g^)+ and is 
accompanied by a term 5(1 — g^) whose coefficient contains -B^a and can be read off directly 
from Eq. fH2]) . 

It is well-known that in the factorization approximation the matrix element B^j^ vanishes, 
and that WA is phenomenologically important only to the extent factorization is actually 
violated. Since Eq. (H2|) is independent of the scale /Xwa, up to 0{as) corrections that we 
neglect, we have in the MS scheme B^j^{fi') = B^^{fi)—p^/2n'^ ln/x'//i. Clearly factorization 
may hold only for a certain value /Xwa = A*/ said, therefore, if factorization holds at /x/, namely 
-Bwa(/^/) = 0, a change of the scale /i/ provides a rough measure of the (minimal) violation 
of factorization induced perturbatively. 

We should stress that the inclusion of WA removes the unphysical singularity in the 
l/mf OPE but not its intrinsic limitations at high g^. The local OPE non-perturbative 
contributions for the rate with arbitrary cuts on Mx, Ei, and g^ is reported in Appendix 
C. It includes WA contributions as discussed above. The main spectra follow from these 
expressions upon differentiation. It can be easily seen that both the rate and some of the 
differential spectra become negative for a sufficiently high cut on g^, a clear sign of the 
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Figure 13: The spectrum in the endpoint region with perturbative corrections switched off. 
The red curve follows from the OPE prediction, the black one modifies it according to Eq. ( fj43l) ) 
to guarantee positivity (6 = 0.4). 



importance of higher dimensional operators. This is apparent in the case of Eq. (HUj) where 
the 1/(1 — (f') singularity drives the spectrum to negative values for q^ > 17.5 GeV^, see 
Fig. [Ill 

Because of the limitations of the OPE in the high region, it is preferable to model the 
tail in a way that is both compatible with the positivity of the spectra and incorporates 

some features of the OPE. We want to modify the OPE expressions before the differential 

rate becomes negative. To this end we change approach for q^ greater than a certain q^ and 

adopt one of the two following methods. 

The first method is based on the local OPE result without the convolution with the 

light-cone distribution functions. One starts from the OPE distribution fHUj) and introduces 

a damping factor in the 1/(1 — q"^) singularity 



dT 
dq^ 



20 f + 66 + 48 r + 74 - 



96(1 



' — hi — 



+ X5(l-g2) + 



(43) 



that maintains the q^ differential rate positive for appropriate values of h (see Fig. [T3l) . We 
have also explicitly written the Dirac delta at the endpoint. The damping factor can also be 
applied to the Ei and Mx differential distributions. In a realistic setting one could expect a 
smooth (positive) bump close to maximal q^. However, since in the following we will always 
integrate over the q^ tail, the rough modelling of (l43l) is sufficient. It should be clear that 
positivity of the rate implies X > 0. The dimensionless parameter X is related to the WA 
matrix element and to the WA scale /i^a by integrating over q^ > In the case of our 



•^This is by far the dominant but not the only negative contribution. Even in the absence of the Darwin 
term the differential rate becomes negative at very high q^. 
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default choice b = 0.4, the minimum value X = corresponds to -Bwa(1 GeV) ~ 0.008 GeV^ 
or equivalently to a factorization scale for WA yU/ ~ 2.2 GeV. The amount of WA implied by 
this ansatz for the high-g^ tail grows with X: X = 0.03 corresponds, in the default setting, 
to SwA(lGeV) ~ 0.017 GeV^ (/i/ ~ 5.7 GeV). Conversely, the experimental determination 
of i?wA(l GeV) would allow us to fix the value of X. 

Since we will mainly be interested in integrated quantities (for instance total rates with 
cuts), this approach is based directly on the formulas in App. C and does not provide a 
triple differential distribution in the high region, but includes the endpoint effect param- 
eterized by X. Perturbative corrections are implemented as for < ql but again there is 
no convolution. We adopt this first method as our default choice. 

In the second method that we employ to model the high region we freeze the distribu- 
tion functions at g^ = g^ and use it in the convolution formula for all > q1- Because of the 
decreasing phase space available at higher g^, this approach effectively amounts to gradually 
phasing all power-corrections out. Although WA does not explicitly enter in this case, the 
ansatz corresponds to a certain amount of WA depending on the value of g^ and on the precise 
values of the other heavy quark parameters. For instance, for g^ = 11 GeV^ and the default 
values of the non-perturbative parameters, it corresponds to i?wA(lGeV) ~ 0.001 GeV'^ 
(/ij ~ 1.1 GeV). An additional amount of WA can be easily accommodated in this frame- 
work as an extra contribution at the g^-endpoint. The triple differential rate here is available 
for each point in the phase space. 

While there is no doubt that the local OPE expansion fails at high g^, the precise deter- 
mination of its range of applicability is an open question. We employ the formulation based 
on the convolution of Eq. fl^ for g^ < ql with 

8.5 GeV^ <ql< 13.5 GeV^ (44) 

keeping the maximum value of g^ lower than the value for which the variance of one of 
the distribution functions gets negative, g^ ~ 15GeV^, because of numerical instabilities 
in dealing with very narrow light-cone functions. To model the high-g^ tail at g^ > q1 we 
employ one of the two methods described above. The difference of the two approaches, as 
well as a variation of the parameters in Eqs. fl42|43ll44p provide us with an estimate of the 
theory uncertainty in the high g^ region. 

6 Results and theoretical errors 

Let us now illustrate our method with a few applications. We start by showing the main 
physical distributions following our default approach, based on Eq. (^11) for g^ < 11 GeV^, 
and on the local OPE modified as in fH3|) for higher g^. We employ the central values of the 
fit in [6j for the non-perturbative parameters at yU = 1 GeV, namely 

mb = 4.613 GeV /i^ = 0.408 GeV^ /x^ = 0.261 GeV^ 
p% = 0.191 GeV^ pIs = -0.195 GeV^ 
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Figure 14: Electron energy (left) and (right) spectra using a few different functional forms. 

1 dr 



r dMx^ 




Figure 15: Invariant hadronic mass spectrum using three different functional forms. 

and set h = 0.4 and X = in Eq. ( l43l) . All numerical results in this section are obtained 
with a C++ implementation of the procedure described in the previous sections. 

Fig. [m shows the electron energy and spectra using a few different functional forms. 
Since we impose the local OPE constraints on the distribution functions at fixed g^, the 
spectrum is by construction independent of the adopted functional form, while the lepton 
energy spectrum shows a limited dependence. The hadronic invariant mass spectrum is 
displayed in Fig. [151 In this case the difference between functional forms is more pronounced. 

One of the checks we have performed concerns the total width that in our scheme is given 
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where Fq = \Vub\'^Gpml{iJ,) /192n^, X = \iamb/2^, and we have left the /x-dependence of mb 
imphcit. The first parenthesis contains the dominant short distance electroweak correction 
[26] which amounts to about 1.4%. The perturbative QCD corrections in (H5|) have been 
known for a while (see [9| for a complete list of references). The non- BLM 0{al) contri- 
butions to the total rate are also known [27] and amount to about +0.5% in our scheme, 
however we do not report them in (1451) as they are absent in our triple differential rate. 
The total semileptonic decay rate can also be obtained by integrating the triple differential 
distribution computed in our default setting with the modified high-g^ tail. In this case the 
result is about 1.4% higher than Eq. fH5l) . using the same input values. This small differ- 
ence is related to the perturbative corrections at the g^-endpoint and might in principle be 
accounted for by a more accurate treatment of the WA contribution, but it will be added 
to theoretical error budget. Using the second option discussed at the end of the previous 
section, namely freezing the distribution functions without taking into account any extra 
effect at the endpoint, the total rate differs from (145|) with default values by —0.8%. 

The next application we consider is the extraction of \Vub\ from some of the latest exper- 
imental results. We consider the following experimental results: 

A Belle analysis with Mx < 1.7 GeV and > 1.0 GeV [2H|; 

B Belle and Babar analyses with Mx < 1.7 GeV, g^ > 8 GeV^ and Ee > 1.0 GeV [281 [29]; 

C Babar with Eg > 2.0 GeV [30] 

This list is far from being complete and it is meant for the purpose of illustration only. The 
analysis B was proposed in [31] and suffers most from the high-g^ uncertainty discussed in 
the previous section (see also Ref. [32]), but in all three cases the high g^ region is probed. 
In comparing with experiment we avoid averaging the various experimental results. In the 
case B, however, we take the arithmetic mean of the two experimental central values as 
reference value. The partial B X^ii^ rate given by the experiments is compared with our 
theoretical predictions, yielding the values for \Vub\ listed in Table 1: 



The central values given in Table 1 refer to our default setting, with an exponential ansatz 
for the distribution functions. The total theoretical errors in the last column of Table 1 are 




(46) 
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cuts 
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exp 


par 


pert 


tail model 


qt 


X 


ff 


tot th 


A [28j 


3.87 


0.71 


6.7 


3.5 


1.7 


1.6 


2.0 


+U.U 
-2.7 


+2.4 
-1.1 




B [281,(29] 


4.44 


0.38 


7.3 


3.5 


2.6 


3.0 


4.0 


+0.0 

-5.0 


+1.4 
-0.5 


±6.61^:1 


C [30J 


4.05 


0.30 


5.7 


4.2 


3.3 


1.8 


0.9 


+U.0 
-6.2 


+1.2 
-0.7 





Table 1: Values of \Vub\ obtained using different experimental results and their experimental and 
theoretical uncertainties (in percentage) due to various sources (see text). / is the estimated 
fraction of events. 



obtained by combining all theory errors in quadrature, with the exception of the asymmetric 
errors due to the functional forms and to X (see below), which are added linearly and kept 
asymmetric. 

We now consider the uncertainty of our theoretical predictions for Table 1. There are both 
parametric uncertainties for instance due to the OPE parameters (m;,, /i^, etc.) and intrinsic 
uncertainties, related to various limitations of our approach. We identify the following 
sources of uncertainty: 

1. the value of a^, m;, and of the other non-perturbative parameters; 

2. higher order perturbative and non-perturbative corrections; 

3. the functional form of the distribution functions; 

4. WA and the high-g^ tail. 

For what concerns the parametric errors (par), we employ ^^(m;,) = 0.22 ±0.02 and take 
the values of the non-perturbative parameters from the fit [6] including all correlations. Not 
surprisingly, the by far dominant parametric error is related to the uncertainty in the h quark 
mass. The value of \Vub\ extracted is quite sensitive to the precise value of the h quark mass. 
For instance, if one employs rrib = 4.677 GeV, as suggested by a fit to charmed semileptonic 
only, instead of rub = 4.613 GeV, the central values in Table 1 change to 3.65, 4.21, 3.79 x 10"^ 
in the A,B,C cases, respectively.. 

To estimate higher order perturbative corrections we i) change the hard cutoff in the 
range 0.7 < ji < 1.3 GeV and ii) rescale the ~ 1% discrepancy in the total rate due to 
perturbative effects in the highest region to the fraction of events. We also consider 
the overall size of the 0(a^/3o) corrections in our analysis (—4.2, —6.0, —6.2% in the A,B,C 
cases, respectively). We take as overall perturbative error in each case the maximum be- 
tween 40% of the 0(a^/3o) corrections and the uncertainty obtained by combining the above 
i) and ii) errors. To estimate higher order non-perturbative corrections and missing per- 
turbative corrections to the Wilson coefficients of power-suppressed operators, we vary the 
non-perturbative parameters within 30% of their central values in an uncorrelated way. This 
leads to errors that are negligible in comparison with those in Table 1. 

The uncertainty due to the functional form (ff) of the distribution functions is estimated 
by comparing a total of about 30 different forms. We display in Table 1 the maximal positive 



24 



and negative variations wrt to the default exponential ansatz. In all cases it amounts to a 
relatively small error. 

Concerning the high-g^ tail we consider three related sources of uncertainty: the mod- 
elling of the tail, the arbitrariness in the choice of the scale where the modelling sets it, 
and the WA matrix element. We estimate the error on modelling the tail by comparing 
our default approach with the second method outlined at the end of the previous section, 
by varying the parameter h in fl43l) and change between 8.5 and 13.5 GeV^. 

In connection to WA we observe that all we presently know about -Bwa is that it is not 
unexpectedly large because no significant excess was found in the high-g^ region by CLEO 
[33]. In the context of our model of the high-g^ tail, this experimental result can actually 
be approximately translated into an upper bound on X, which in turn is related to i?wA oi 
11 f. The CLEO bound is not very strong, X < 0.07 at 90%CL, and still falls short of the 
expected size for the WA contribution to the total rate (5rwA/rtot ~ 2% [3l]. We vary X 
in the conservative range < X < 0.04. In terms of model-independent parameters this 
corresponds to 0.008 GeV^ < 5wa(1 GeV) < 0.020 GeVl In the future, it will be possible to 
measure the size of the WA expectation values using, for instance, measurements of the g^ 
spectrum close to the endpoint and of its moments [23, and by comparing the semileptonic 
rates of and . 

In our default set-up we have assumed X = 0, but while we expect it to be positive and 
possibly small, there is no compelling reason for it to vanish. The X error in Table 1 should 
therefore not be treated as gaussian. It is worth observing that a non-zero value for X would 
bring the different experimental determinations of \Vuh\ closer to each other. For instance, 
X = 0.04 would lead to iKbl = 3.76, 4.21, 3.80 x 10"^ in the A,B,C cases, respectively. 
Suggestive as this may sound, it is clear from our analysis that an upper experimental cut 
on g^ would improve the theoretical precision in the inclusive determination of \Vub\- In the 
absence of a dedicated experimental analysis, we can subtract Belle's results for the cases A 
and B [28] and extract a value of \Vub\ for the combined cut Mx < 1.7 GeV, Ei > 1.0 GeV, 
and g^ < 8GeV^, i.e. with an upper cut on g^. The cut on g^ is relatively low and the 
predicted fraction of events in this case is only 33%, but the measured rate is lower than 
expected: it corresponds to \Vub\ = 3.18 x 10^'^ with -^A.b'^^Vo theoretical error, dominated 
by parametric uncertainties. We were unable to compute the experimental uncertainty for 
the g^ < 8 GeV^ case, which will certainly be larger than those in Table 1. Such a low 
value of \Vuh\ may, however, indicate either an experimental problem or an underestimate of 
theoretical errors in the high g^ region, which is common to the three cases considered in 
Table 1. 

7 Conclusions 

We have calculated the triple differential distribution for inclusive semileptonic decays with- 
out charm, B X„£z/, consistently including all the known perturbative and non-perturbative 
effects, through 0(a^/5o) and 0(l/m^), respectively. Our theoretical framework is based on 
the OPE and incorporates a hard Wilson cutoff ~ 1 GeV. This involved new perturbative 
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calculations discussed in Section 2. 

Our approach has several new elements that we have listed in the Introduction and 
explained in detail throughout the paper. We recall the main ones: we parameterize the 
Fermi motion in terms of a single light-cone function for each structure function and for 
any value of g^, and include consistently the sub leading effects; this is accomplished at the 
same level of model-dependence as for the leading twist distributions; we implement for the 
first time the complete BLM corrections to the triple differential rate; we present a detailed 
discussion of the high-g^ tail and of Weak Annihilation effects. Our approach is completely 
implemented in a C++ code that is available from the authors. 

We have extracted \Vub\ from some of the latest and most precise experimental data, 
providing a detailed estimate of the theoretical uncertainty. Our results, hsted in Table 1, 
agree within theoretical errors with those obtained with other methods used by HFAG [121 
[T3] . We find that the dependence on the functional form assumed for the distribution 
functions is rather weak. However, the critical role played by the high-g^ tail becomes 
evident from our analysis; it contributes in a significant way to the theoretical uncertainty 
of the present inclusive determinations of \Vub\- We have modeled the high-g^ region in two 
different ways complying with the positivity of the differential spectra and have accounted 
for the WA contributions. We find that non-vanishing WA effects tend to suppress the value 
of \Vub\ extracted from the data and, using recent Belle results, we have argued that the 
low-g^ sample of events leads to a lower value of \Vub\ that conflicts with that extracted from 
the g^ > 8 GeV^ sample. 

Since the high-g^ tail presently leads to a sizeable uncertainty, we encourage our exper- 
imental colleagues to pursue analyses with an upper cut on g^, and to perform an accurate 
separate measurement of this domain. 
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Appendices 

A Perturbative corrections 

In this Appendix we provide all necessary analytic expressions for the perturbative correc- 
tions to the three structure functions in the scheme with a hard cutoff /i. In the absence of 
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cutoff the 0{as) real gluon emission terms can be gleaned from [T5| [9]: 
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where u = 1 — 2go + g^ and 
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The plus distributions in Eqs. fj471 - l49|) are defined by (n > 0) 
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where G(go, g^) is a smooth function. 

The NLO real emission contributions to the structure functions in the presence of 
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Wilsonian cutoff, see ([7]), are a new result and read as follows: 
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We also report the soft- virtual structure functions in the absence of the cutoff [15], [9] . 
They are consistent with the way we have performed the subtraction in the real emission 
contributions. At 0{as) they are 
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while at 0(a^/3o) they are given by 
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B Structure functions in the local OPE and their 
go-moments 

In the adopted normalization, the power corrections to the structure functions read 
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We have used the short-hand notation 5„ = 5*^")(1 + g^ — 2go), where the n-th derivative of 
the Dirac delta is taken wrt its argument. We give explicit expressions for the zeroth, first 
and second go-moments, at fixed g^, of the three form factors, up to C(l/m^) corrections in 
the OPE. Separating the tree-level and power corrections contributions, we define 
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and 
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The functions /'s defined in Sec. |3] correspond to central moments and are linear combi- 
nations of above M^-''^'^^^'^ or M-*-"''*'^""' moments: 
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The explicit expressions for the zeroth moments are: 
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The explicit expressions for the second moments 
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C Local OPE results with arbitrary cuts 

In this Appendix we report the local OPE results for the rate of S — > X^lv subject to 
standard cuts on Mx, E£, and q^. The expressions contain only the non-perturbative power 
corrections. We adopt the following notation 
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where E^^cut and q^^^ are lower cuts on E^^^q^, and M'j^^^^^ is an upper cut on the invariant 
mass of the hadronic system. Only r, the effective lower cut on g^, is relevant for power 
corrections: the upper cut on the hadronic mass must be lower than VAMb ^ 1.8 GeV to 
play a role in the following expressions. There is also a relation between ^ and r: both cuts 
are relevant only when r < ^ < 1, otherwise only r is relevant. The total width subject to 
the above cut is given in the two regions by 
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and 
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